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1 Introduction 



Since the historical experiments of the SLAC-MIT group [jT| revealed the partonic character of 
the nucleon's constituents, deep-inelastic scattering (DIS) has been one of the most important 
tools to increase our understanding of the inner structure of the nucleon. 

The first experiments established that nucleons can be viewed as consisting of point-like 
constituents, which later were identified as quarks and gluons. The observed scaling behavior 
of the lepton-nucleon scattering cross section was interpreted as incoherent scattering of the 
probing lepton with the partonic constituents of the nucleon. 

Since the advent of quantum chromodynamics (QCD) considerable substance has been added 
to this simple picture. The scaling behavior has found an explanation through the asymptotic 
freedom of QCD. Asymptotic freedom combined with the factorization property resulted in a 
systematic expansion of the DIS cross section in terms of the coupling constant as{Q'^), evaluated 
at the scale set by the virtuality of the hard photon. The leading order term in this expansion 
accounts for the simple parton model result, while higher order corrections modify the scaling 
behavior, resulting in a logarithmical scale dependence. 

The QCD explanation of logarithmic scaling violations in DIS, as observed in a later set of 
experiments 0, is considered as one of the most important successes of the theory. 

Since then, the experimental precision of DIS measurements has reached remarkable accuracy 
over a wide range in Q^. Using these data it has been in principle possible to study not only the 
logarithmic corrections to scaling, but also to extract higher twist terms, i.e., power suppressed 
corrections that fall off like powers oi . 

From the theoretical side, the best tool to analyze such power corrections is the framework 
of the operator product expansion (OPE) Twist-4 corrections to DIS have been studied 
systematically e.g. in |^. Kinematical corrections arising from the non- vanishing mass M of 
the target hadron fall off like powers of M'^/Q'^ and can be attributed to power suppressed 
contributions of twist-2 operators. They can be taken into account exactly by introducing the 
so-called Nachtmann scaling variable ||^. 

The other set arises due to the higher twist operators that are sensitive to multi-parton 
correlations in the target. While estimates of twist-2 matrix elements from lattice QCD are 
already available 0, reliable estimates of higher twist matrix elements are not yet feasible. 

In particular, calculation of higher twist operators with quantum numbers which do not 
prohibit mixing with lower twist operators is a severe theoretical problem that has not yet been 
solved [g. 

The problem originates from the fact that twist-4 operators, in addition to the usual loga- 
rithmic scale dependence due to their renormalization, may exhibit quadratic UV divergences. 
Recall that the matrix element of a twist-2 operator is a dimensionless number. A related twist-4 
operator of the same spin and quantum numbers has therefore a matrix element of dimension 
2. Then, radiative corrections result in a contribution of the form of a square of the UV cut-off 
multiplied by the lower order matrix element of the twist-2 operator. Such a mixing makes the 
definition of the twist-4 contribution ambiguous. In the OPE of a physical quantity, like DIS 
structure functions, this ambiguity always cancels against the corresponding ambiguity in the 
definition of the twist-2 contribution. The latter arises because of the asymptotic character of 
the QCD perturbation series |10]. Hence, the sum of twist-2 and twist-4 contributions is un- 
ambiguous up to order l/Q^, provided that both are calculated within the same regularization 
scheme. 

Recently, this subtle relation between twist-2 and twist-4 contributions has motivated a 
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phenomenological hypothesis [T^ stating that the main contributions to matrix elements of twist- 
4 operators are proportional to their quadratically divergent parts^. In processes which cannot 
be analysed in terms of OPE one extracts information about power corrections directly from 
the large-order behavior of the corresponding perturbative series. Power suppressed corrections 
to various observables, like event shape variables in DIS [|12|] and e~ annihilation, as well as 
DIS structure functions [|13], 0, [l^ have been shown to follow the behavior predicted by the 



UV-dominance hypothesis (for recent reviews see ||T6|, [T^])- Hence, with all reservations, it can 
be considered as a useful phenomenological tool for estimating power suppressed contributions. 

Comparison with the existing experimental data on leading power corrections to the F2{x, Q^) 
structure function has shown that the non-singlet IR renormalon calculation describes the 
data on proton and deuteron structure functions very well in the region of Bjorken-x > 0.25 but 



there is a systematic discrepancy in the region x < 0.25 [T^. In the present paper we extend the 
existing analysis of power corrections to flavour non-singlet structure functions |]TB|, |T^] to 
the quark pure flavour singlet case. We follow the idea to trace twist-4 operators by calculating 
their UV divergent part in the renormalon approach. As expected, the resulting correction turns 
out to be much smaller than the non-singlet one in the large-a: domain, but it is substantial in 
the small-x region. It is encouraging to see that our predictions agree with the tendency seen in 
the data points below x < 0.2 0. 

Extending the existing measurements deeper into the small- a; - small-Q^ domain of DIS one 
ultimately enters a transition region between interactions of hard and (almost) real photons with 
a nucleon. The twist expansion provides a tool to approach this transition from the large-Q^ 
side. The importance of twist-4 corrections to F2{x, Q"^) in the small x domain has been realized 
a long time ago [|1^ . The four-gluon operator was identifled as a potential source of large twist-4 
corrections [|19], and its leading log anomalous dimension was calculated |]19|, |2^. Recently, a 
model analysis of twist-4 contributions has been performed in the diagrammatic language 
We emphasize that the renormalon analysis of twist-4 contributions in the pure flavour singlet 
sector, presented in the present paper, cannot substitute a non-perturbative QCD calculation of 
these corrections. Nevertheless, as we shall discuss in details below, certain interesting features 
arise which can be confronted with existing data. 

Recently, new data have been published ^ which extends a previous NMC analysis 
[El, EH] further into the small-x - small-Q^ region. Down to of the order of 1 GeV^ the data 



can be well described by a set of radiatively generated twist-2 parton distributions ||26[. On 
the other hand, even if the precise form of its x-dependence is still subject to debate, there is 
certainly no reason to assume that twist-4 corrections are small and that they can be neglected. 
In our opinion, it would be interesting to have the same data reanalyzed using a model for the 



twist-4 contribution e.g., in the form derived in the present paper or taken from Ref. [21 



In the following we will investigate power corrections to the three structure functions F2, Fl 
and Qi that appear in the well known decomposition of the hadronic scattering tensor of deep 
inelastic lepton nucleon scattering p7| : 



W,,{p,q) = ^ [d'ze''^{pS\Uz)UO)\pS) 

ZTT J 



<hi<iu 

q2 



—Fiix, Q^) - {gf,u + P^Pu 



2x 



{p-qf 



p-q 



2x 



^ If quantum numbers prohibit mixing of the twist-4 operator with a lower dimensional twist-2 operator, the 
former exhibits, of course, no quadratic divergences of its matrix elements. It would be very interesting to find 
a set of experimentally accessible observables which would allow to extract a power correction which cannot be 
interpreted as a UV dominated twist-4 matrix element. 
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-ief^uXa gi{x,Q'^) . (1) 

p-q 

Here is the electromagnetic quark current, x = Q'^/{2p ■ q) and q^ = —Q^. The nucleon 
state \pS) has momentum p and spin S with S*^ = — M^, M being the nucleon mass. The 
convention for the e-tensor has been taken from [^. We have neglected terms arising due to 
weak interactions as well as the second spin dependent structure function g2{x,Q'^) which is 
suppressed kinematically. 

Our presentation is organized as follows. In section 2 we present the basic assumptions and 
definitions of the renormalon model. In section 3 we describe our calculation, which leads to 
the estimate of the x-dependence of the twist-4 corrections to DIS in the flavour singlet channel. 
Explicit calculation is performed for the pure-singlet quark contribution, and subsequently this 
result is used to model the x-dependence of the corresponding gluon contribution. Section 4 
is devoted to phenomenological discussion of our results. We discuss the renormalon model 
predictions for 1/Q^ corrections to unpolarized and polarized nucleon structure functions with 
a particular emphasis on the small x-dependence of these corrections. A possible data fitting 
procedure is discussed in the summary. Finally, our appendix contains the list of all renormalon 
model formulae derived in the present paper. 



2 Nucleon structure functions beyond the leading twist 

According to the OPE, hadronic structure functions Fj, i = L, 2 in can be decomposed up to 
O(^) accuracy as 

F,{x, Q') = Fr-'{x, Q') + Q') + ^,K{x, Q') + O(^) , (2) 

where F^'"~'^ describes the leading twist-2 contribution, h™'^ describes the target mass cor- 
rections which are directly related to twist-2 matrix elements ||2^. It is hi which contains the 
genuine twist-4 contribution and is sensitive to multi-parton correlations within the hadron. 
The goal of this paper is to provide a phenomenological model for coefficients hi{x,Q'^) in the 
flavour- singlet sector. 

In the QCD-improved parton model the twist-2 contribution to a deep inelastic structure 
function can be represented as a convolution of process independent, universal parton densities 
with perturbative coefficient functions. For a general number of flavours Uf the corresponding 
formula reads: 



r-'ix,Q') = X - (-Eel 



1 d^ r 1 



+G{x/ z, f?) Ci^ciz, Q V 



(3) 



Here G(x,yU^) denotes the gluon density, S(x, /i^), A(x, /i^) stand for the singlet (S) and non- 
singlet (NS) combinations of quark densities, and Ci^c, Clg, and Cj^^^ represent Wilson coeffi- 
cients in the corresponding channels. represents Fl or F2. The factorization scale is denoted 
by /i^. 
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The flavour singlet combination of the quark densities is defined as 



(4) 



where fi and fi stand for quark and antiquark densities of species i. The non-singlet combination 
is given by 



el- 



i=l 



{Mz,f,') + Mz,fi')) 



(5) 



The charge of the quarks is denoted by e,. A decomposition completely equivalent to Eq. (|^) 
can be written for the polarized structure function gi with polarized quark and gluon densities 
AS and AG, and coefficient functions ACg and ACg, respectively. 

The parton densities are pure twist-2 quantities. Taking moments of structure functions one 
obtains the OPE relation 




Here a<; stands for 



9' 



167r2 



(6) 



(7) 



and the An are the matrix elements of the spin-N twist-2 operators which determine the parton 
densities in Eq. (H). Fk can be Fl, F2 or x^^i.The coefficient functions can be calculated order 
by order in perturbation theory. At present, Wilson coefficients are available up to the second 
and in some cases, for the lowest moments, up to the third order in ||3l[ . 



order 30 



Matrix elements of twist-4 operators |^ contribute power corrections to the simple partonic 
picture in Eq. (|^). According to the UV dominance hypothesis they can be considered propor- 
tional to their quadratically divergent part or, equivalently, to the uncertainty of the perturbation 
series in the definition of the twist-2 contribution. Commonly, the divergent series of radiative 
corrections is regarded as an asymptotic series and defined by its Borel integral. The actual 
calculations are done in the limit of large nj, which allows to resum the fermion bubble-chain 
to all orders yielding the coefficient of the a^n^~^ - term exactly. Subsequently, it is converted 
into the exact coefficient of the a^PJ^~^ - 
known as the 'Naive Non-Abelianization' (NNA) 



term by the substitution Uf ^ Uf — 33/2 = — |/?o; 



2], The asymptotic character of the re- 



sulting perturbative series leads to resummation ambiguities - a singularity in the Borel integral 
destroys the unambiguous reconstruction of the series and shows up as a factorial divergence 
of the coefficients of the perturbative expansion. The general uncertainty in the perturbative 
prediction can be estimated to be of the order of the minimal term in the expansion, or by taking 
the imaginary part (divided by tt) of the Borel integral. Both procedures lead to resummation 
ambiguities of order (A^/Q^)*^, with r = 1 for the leading IR renormalon. The resulting model 
for the twist-4 contributions hi, see Eq. (^, has the form of a Mellin convolution : 

KG{z,^^') = {A.,,Yx [' -A^l{z,\og{{A.,GY/f^'))G{x/z,f,') 

J X Z 

hU^,f^') = {Atrx [' -Alg^'\z,\og{{Afjyf,'))J:{x/z,f,') 

Jx Z 
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K^{z,^^') = (Af/rx/'^</(^)(z,log((Af/)V/i^))A(xM^^), (8) 

J X Z 

with G, S and A being the twist-2 parton densities of Eq. (^. The scale dependence of the 
coefficients hi is at variance with the QCD predictions, as it is related to anomalous dimensions 
of twist-2 rather then twist-4 operators. In order to escape this difficulty we suggest that Eq. (§) 
should be understood as a phenomenological model valid at some low scale [j? of the order of 
1 — 2 GeV^, where higher twist contributions become important. The mass scales Aj ^ have to be 
fitted to experimental data. The non-singlet coefficients ^j^'^'' were calculated in []T^ |T^, n, 

The pure-singlet coefficients J^f* (and corresponding coefficients J^i^ ^ related to or twist-6 

corrections), which are the main result of the present paper, have been collected in the appendix. 



3 Calculation 

In this section we describe the calculation of the coefficients J^^^ and J^^^ ■ They can be obtained 
in many ways, either by calculating one-loop Feynman diagrams with an infrared regulator such 
as a gluon mass or using a dispersion representation of the 'effective' coupling constant 



discussed in [|T2[, or by analysing the large-order behavior of the twist-2 coefficient functions. 
We follow the third method by investigating the large-order behavior of the coefficient functions 
in the limit that the number of flavours, i.e. rij, goes to inflnity. In the case of modeling power 
corrections to NS-structure functions it is sufficient to compute one-loop diagrams, while the 
evaluation of the singlet part is more complicated and involves two-loop calculation. The ef- 
fective coupling method has been recently used in the two-loop calculation of power corrections 
to photon-photon scattering by Hautmann Up till now, the equivalence of all three afore- 



mentioned methods of computing A!f^ and has been firmly established only at a one loop 
level. In case of a two-loop calculation like the present one the equivalence between renormalon 
approach and the flnite gluon mass method persists, as discussed below, but the relation to the 
effective coupling approach |jl2| has not yet been proven. 



The singlet coefficient may be decomposed as 

where PS denotes the so-called 'pure-singlet' part where the photon line is attached to a closed 
quark loop which is not sensitive to the flavour of the target. While the large- riy expansion of 
Cf^ is already known, here we describe the large-rij expansion of , see diagrams in Fig. ^ 
To investigate the large-n / behavior of the gluonic coefficient function Ci^c one has to analyze 
fermion box diagrams with external gluon legs where in addition a fermion bubble chain is 
inserted, see Fig. ^. Explicit enumeration of the rif factors arising in both cases shows that 
contributions of the latter type are suppressed by a factor of 1/n/ as compared to the quark 
pure-singlet part. The large-n/ structure of the quark pure-singlet coefficient can be written as 



i,q 



(10) 



The curly brackets indicate the contribution in the nj expansion that comes from the fermion 
box. The gluonic coefficient is formally of the same order in nj. Taking into account that the 
lower vertex couples to the quark singlet or gluon structure function one obtains effectively an ad- 
ditional factor of Uf in the former case. The quark singlet structure function (^ J2kLi efc)S(x, /i^) 

is 0{nf) as compared with the gluon structure function G^x.jj?)^. Recall, for example, that in 

^We thank V.Braun for useful discussions about this point. 
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Figure 1: Diagrams which contribute to the 'pure-singlet' part of the quark coefficient function. 



the large-n/ limit the fraction of longitudinal momentum carried by quarks and gluons is of the 
order 1 — 0{l/nf) and 0{l/nf), respectively |^6|]. Following this argument one has to calculate 



only the quark singlet part as the dominating contribution in the large-n j limit. 

Instead of calculating the corrections term by term to obtain the large-order behavior of 
the perturbative expansion, it is convenient to deal directly with the Borel transform of the whole 
series. In the special case we are considering here, there are two bubble chains which contribute, 
both having the same virtuality — fc^, see Fig. |l|. By performing the substitution a""*"^ u^/nl 
one easily obtains the Borel transform of the square of the running coupling 



BT 



,l+/5oa,log(-fcV(/i2e-^)), 



with s = PoU, and u being the Borel parameter. The Borel transform can be viewed as generating 
function for fixed order coefficients. Since our expansion starts at order the Borel transform 
vanishes at s = as it should be. Inserting Eq. (|ll]) in the loop calculation this simple picture can 
be destroyed when collinear poles manifest itself as 1/s contributions which effectively simulate 
finite contributions to the Borel transform. Those poles are first to be subtracted to allow for 
the determination of fixed order results. 

At this point it is convenient to discuss the equivalence between the renormalon approach 



and the finite gluon-mass method, where functions Af^ and A\'^' are identified as coefficients of 
non-analytic terms in the small gluon-mass expansion of diagrams in Fig. |^. Note that in the 
present case both gluon chains depend on the same virtuality k"^. Following ||3^, one can derive 
the Mellin representation of the square of the massive gluon propagator in the form: 



l(4) 



1 



1 1 



-1/2+ioo 



(A;2 - A2)2 27ri A;^ 



du 



7r{u + l) ( \^ 



Sm VTM 



(12) 



-l/2-ioo 
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Figure 2: Typical diagram which contributes to the gluon coefficient function. 



which can be rewritten as 



-l/2+joo 



(A;2 - A2)2 27ri 



du 



l/2-ioo 



Msin(7rM) \/i^e' 



-c 



BT[as{k% 



(13) 



see Eq. (|TT|). Inserting this equation into Eq. (|T6]) below and inverting the resulting formula 
by means of the Laplace transformation, it is easy to see that the one-to-one correspondence 
between renormalon poles and non-analytic terms in the small gluon-mass expansion, discussed 
holds at the two-loop level as well. On the other hand, relation between the renormalon 
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m 

method and the effective coupling approach still requires further investigation [^. 

To obtain the IR-renormalon poles from the diagrams in Fig. |I| we found it advantageous 
to follow the procedure for calculating DIS coefficient functions explained in ||30|, 0. In this 
method the hadronic scattering tensor is calculated directly as the cross section of photon-parton 
scattering. Recall that in the partonic picture the hadronic structure functions can be written 
as 

"1 dz, 



i^*-2(a;,g2)^^^ / -ffc(x/z)^,fc(z,Q^£) 
1. Jx Z 



(14) 



where f^ and jFj ^ denote the bare parton density and the parton structure function corresponding 
to a parton k = q,q, G. The parton structure function can be obtained order by order in the 
perturbation theory by calculating the radiative corrections to the parton subprocess 



l\q)+k{p)^X 



(15) 



where a virtual photon of momentum q hits a parton k with momentum p. In the general case, 
Ti^k contains collinear divergences. They can be regularized e.g., by performing the calculation in 
d = 4: — 2e dimensions. The dependence of !Fi^k on the regulator is indicated by the parameter e. 
The infrared-divergent part has to be subtracted from the partonic cross-section and factored into 
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the bare parton density to define a physical, measurable quantity. In the present calculation it is 
the Borel parameter s which alters the power of the gluon propagator, and therefore regularizes 
divergent integrals. Because of that, it is sufficient to perform the calculation in ci = 4 dimensions. 
The collinear 1/e poles then show up as 1/s poles. Since we are only interested in the power 
corrections which manifest themselves as poles at s = 1,2,..., there is no need to calculate 
the subtractions explicitly. Note that the anomalous dimensions in the MS scheme do not 
exhibit renormalon singularities, and therefore the subtraction terms do not introduce additional 
renormalon poles |38|, |4l|. Hence, the IR-poles in the Borel image of structure functions Fi{x, Q^) 
arise only from the IR renormalon poles in the Borel image of partonic structure functions J-'i^k, 
calculated without subtractions. 

Partonic structure functions can be obtained from the partonic tensor 

1 1 oo „ 

W,Ap, = 2 E y dPS^'^M,{l)M,iiy , (16) 

where / dPS^''^ represents the /-body phase space integration. The factor ^ comes from aver- 
aging over quark spins in the unpolarized case. The tensor Wfj_iy{p,q) allows for an equivalent 
decomposition in terms of J^2,k, ^L,k, gi,k as the hadronic scattering tensor in Eq. (0). Choosing 
appropriate projections it holds: 



12f 

1 



= - m - -^p.p.w^'' 



gi,k = e.^x.q'p'^W^'' . (17) 

2 q ■ p 

The normalization is chosen such that in the simple parton model we get 

^iyz,Q^) = 5(1 - ^) ,^i^J(^,Q2) = a,CF{Az) , (18) 
3.1 Quark coefficient functions 

To compute the pure-singlet contribution one needs to calculate the diagrams in Fig. |l[ In the 
improved parton model this corresponds to the calculation of the two-to-three-body subprocess 

l*{q)+p-^h + k2 + p' , (19) 

which requires the three body phase space integrals 

dP5(3) = -J— / d''kid''k2d''p'5^^\kl)5^+\kl)5^+\p''') (20) 
(27r)"' J 

x6{p + q - ki - k2 - p') 

that appear in Eq. (pIBf) . It turns out to be convenient to evaluate the three particle phase space 



in the CMS system of the outgoing qq pair of the fermion box. Introducing variables as in 
we find for the phase space integral 

1 „ol -2 . . 



dPS^'^ = -—Q' -f dy{l-y)( dv r d<P r de sm{9) . (21) 

[AlT)^ z Jo JO JO JO 
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Here 9 and are the polar and azimuthal angles in the CMS system of fci and k2 respectively, 
while V and y are rescaled u and t channel invariants. The exact definition of theses variables 
can be found in the appendix of For our purpose the definition of y 



y 



■ z 



(22) 



is important, where —k^ is the virtuality of the dressed gluon and z is the argument of the partonic 
structure function j^ij.. It follows that the Borel transforms of dimensionless projections of the 
partonic tensor Eq. (|17D can be written in the form 



Jo z \y fio 



^{z,y) 



(23) 



which shows that the dependence on the Borel parameter s factorises and is infiuenced only by 
the y integration. To find the renormalon poles it is therefore not necessary to evaluate the whole 
expression exactly. Instead, we can expand J^{z, y) in y. The integration over y transforms terms 
~ y"' into simple poles ~ s/(?t, — 1 — s) and terms ~ y^\og{y) into double poles ~ s/(?t, — 1 — s)^. 
Note that \og{y) terms arise from the integration over the polar angle 9 where partons ki and 
/c2 become coUinear. 

The expansion of T{z^ y) in y then leads to an expression of the form 



:F{z,y) = Y: iAr,iz)+B^iz)\ogiy)) 



n=0 



(24) 



The expansion coefficients An{z) and Bn{z) can be expressed as linear combinations 



An{z), Bn{z) ~ \0g{z)Wi^n{z) + Wo^n{z) 



with the coefficients Wi^n{z) and Wo,n(^) being simple polynomials in z. Inserting Eq. (p^ back 
into Eq. ( |23| ) we observe that although the form of the phase space integral suggests that the 
Borel transform BT[J-'{z)]{s) is simply proportional to z^, after performing the y integration one 
actually obtains an expression of the form 



yS + 1 



const. + (Aoiz) - Bo{z)) ^ Bi{z) — 

1 + s s 



+ iA2iz) + B2iz)) 



z 



'1 



yS-1 



1 - s 

,s-2 



+ {2A,{z)+B,{z)) 



z-' 



2-s 



Bo(z) 



2B.(z) 



{2-s) 



O 



3-s 



(25) 



It is interesting to note that the coefficient of the 1/(1 + s)-pole in the sum over all diagrams 
vanishes, i.e. that the coefficient function of a DIS structure function does not contain any UV- 
renormalon poles. This is in line with previous observations that positions of UV-renormalon 



poles in NS-structure functions always depend on the Mellin moment |T3[ , and that therefore 
UV-renormalon contributions are absent after the inverse Mellin transformation is performed. 

Note the appearance of double poles in Eq. (^), which were not present in one-loop calcula- 
tions of renormalon contributions to DIS structure functions in the non-singlet sector [0, 0, 15 . 
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By considering renormalisation group improved version of OPE it can be shown that renormalon 
singularity of the Borel transform is in general expected to have the form 



BT[Hz)]{s) ~ , ''^^f^- . (26) 

^ ^ '^^ ' (n — s)l+7o//3o ^ ' 

Here 70 corresponds to an eigenvalue of one-loop anomalous dimension matrix of higher twist 
operators of dimension n which contribute to OPE of the hadronic tensor in Eq. (Q). Calculation 
of anomalous dimensions considerably simplifies in the large-n/ limit, resulting in eigenvalues 
which are either equal to zero or to a multiple of (3q. The double pole found in Eq. (p5|) arises 
from operators with 70 = /?o- 

Terms in Eq. ([25| ) singular at s = are directly related to singularities which occur in a 
calculation of J^{z) in = 4 — 2e dimensions. The coefficient of the 1/s pole corresponds to the 
quadratic collinear divergence and is directly related to the pole in the usual dimensional 
regularization. In general, this coefficient can be written as a convolution of the leading order 
splitting functions 

BAz) PjS^ ® P!S> . (27) 
Thus, for the structure function F2 one obtains 

Bi{z) ~ Pg(J) ® P^(J) = CpTpUf^ [4(1 - xf + 3x(l -x) + 6(1 + x) log(x)] , (28) 

while for gi one finds 

Bi{z) ~ AP^^J) ® AP^(J) = 32 CpTpUf [5(1 - x) + 2(1 + x) log(x)] . (29) 

Note that for the longitudinal projections the fermion box is free of collinear divergences and 
therefore the contribution does not contain any pole. 

The form of the constant term in the expansion can also be determined from general consider- 
ations. It is given by a combination of next to leading order splitting function and a convolution 
of leading order coefficient function and leading order splitting function corresponding to single 
1/e poles. However, determination of the constant term in our approach requires the expansion 
to all orders given in Eq. (p^). 

Recall that for x close to 1 one expects the ratio of twist-2 to twist-4 contributions to behave 
like h?/{Q'^{l — x)), or approximately like A? /W'^ The same trend at small x would 



correspond to an expansion in terms of x / Q^, see e.g. Ref. . However, the exact calculation 
which takes into account Eq. ( P^D shows that in the small x region the expansion parameter, see 
Eq. (|31|) below, turns out to be A^/Q^ instead of x /Q'^. As a consequence, in the renormalon 
model the ratio of the twist-4 and twist-2 or the twist-6 and twist-4 contributions goes to a 
constant limit as x ^ 0. Note that in the renormalon model calculation of power corrections to 
parton fragmentation functions in e"*" e~ anihilation into hadrons the expansion parameter turns 
out to be A^/(x^Q^). Hence, the expansion breaks down in the region x < A/Q and has to be 
resummed [^. No such resummation is required in the present case. 

To obtain the coefficients Ai^k in Eq. (^ one has to transform the Borel images back to the 
as representation by using the inverse Borel transformation 

1 °° 

:Fiz;as) = — I cise-^/(^°"^)PT[^(z;a,)](s) . (30) 
Po i 
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An unambigious reconstruction of the perturbative expansion is prevented by the poles on the 
integration contour and therefore the summed series acquires an imaginary part which is directly 
related to the UV-divergent part of the higher twist operators we are looking for. Its magnitude 
can be estimated either directly or by taking the imaginary part (divided by vr) Q/n of the 
inverse Borel transform in Eq. ([30| ) P6|| . This procedure leads, for the simple and quadratic 
poles of the right-hand side of Eq. ( p5|) , to 

The ambiguity in the sign of the IR-renormalon contributions is due to the two possible contour 
deformations, namely above or below the pole aX s = n. 



3.2 Gluon coefficient functions 

As it has been mentioned above, the gluonic contribution is subleading in the large-nj limit. For 
finite nj, however, gluons certainly play an important role. To model their contribution we have 
adopted a procedure introduced in in the analysis of power corrections in fragmentation 
processes in e^e~ annihilation. The idea is to represent the pure-singlet unpolarized quark 
coefficient function Af^ as a convolution of the leading order quark-gluon unpolarized splitting 
function 

p(;)(.) = c,l±il^ (32) 

with an effective gluonic unpolarized coefficient function Ai^ciz)- After taking Mellin moments 
of both sides of the equation 

A^a ® = A^f , (33) 
Aifi can be obtained by performing an inverse Mellin transformation of 

AS (AT) 

The same procedure can be applied also to define the polarized gluon coefficient function Aa,g{z), 
with the obvious modification that now the polarized pure-singlet quark coefficient function and 
the polarized quark-gluon splitting function 

AP^'J{z)=4Cp{2-z) (35) 

enter Eq. (p4|). The obtained formulae are presented in the appendix. In Fig. ^ and Fig. ^ we 
have compared the ratios of quark and gluonic twist-4 corrections to F2{x, Q"^) and to gi{x, Q^), 
respectively. It turns out that in both cases the shape of the gluon contribution which results 
from deconvolution procedure described above is quite similar to the shape of the original quark 
contribution. It is also easy to show that the ratio of quark and gluon contributions tends to a 
constant if a; — 0. As the overall normalisation of each contribution is anyhow an adjustable fit 
parameter, we have modified the ratios in Fig. ^ and Fig. |^ by dividing in each case the quark 
contribution by a factor 2. 
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Figure 3: Ratio of the quark pure-singlet contribution to the gluonic contribution obtained according 
to Eq. ( |34D for F2. The dotted line shows contributions of the single poles alone, the solid line 
corresponds to the coherent sum of contributions from single and double poles. We have adjusted 
the normalization by dividing the pure-singlet quark contribution by 2. 



Note that this procedure results in the twist-4 contribution originated by gluons which de- 
pends hnearly on the twist-2 gluon longitudinal momentum distribution. On the other hand, it 
has been known for a long time that an important class of twist-4 corrections to DIS originates 



from matrix elements of gluonic operators involving four gluon field strength tensors |T^, thus 
being proportional to the square of the gluon momentum distribution. This is a good illustration 
of the fact that the renormalon model does not provide a description of all contributions to power 
suppressed corrections, and that some of the missed terms can be important for phenomenology. 



4 Discussion 

As it is known that the renormalon model cannot account for the absolute normalization of 
the twist-4 corrections, we follow the attitude proposed in 



11, 18 , and assume that model 



parameters, the overall normalization factors and mass scales, should be fitted to the data. 
Recall that in the flavour non-singlet case data at x > 0.25 can be fitted using the normalization 
factor about 2 times larger than = -^A^e"*" arising directly from the renormalon model. 
To get a feeling about the shape of pure flavour- singlet twist-4 corrections, in the following we 
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Figure 4: Ratio of the quark pure-singlet contribution to the gluonic contribution obtained according 
to Eq. (15) for gi. The solid line shows a coherent sum of contributions from single and double 
poles. We have adjusted the normalization by dividing the pure-singlet quark contribution by 2. 



discuss the renormalon model predictions for twist-4 corrections to F2{x,Q'^), gi{x,Q'^), and 
for a deuteron target, obtained at a fixed mass scale. For definiteness, we have used 
the NLO MRSA parameterization |^ and the NLO Gehrman-Stirling gluon-A set p8[ at = 4 
GeV^ as an input for twist-2 unpolarized and polarized parton distributions, respectively. We 
have checked that other sets of LO or NLO parton distributions produce similar results. While 
such an analysis cannot substitute for a data fit, it transparently reveals the gross features of 
the model. 

Fig. I shows the prediction for the ratio Fl'^-^{x,Q'^) / Fl^-^{x,Q'^) at = 4 GeV^ for a 
deuteron target as compared with the data from Ref. To obtain this plot we have fixed the 
sign in the right-hand side of Eq. (0), such that contributions arising from both simple and 
quadratic poles add coherently. At large x the renormalon model prediction is dominated by 
contributions from NS-type graphs, common for singlet and non-singlet channels, see Eq. (^). 
It explains why it has been possible to fit deuteron and proton data in this domain by the 
renormalon model correction based only on the non-singlet coefficient function \Ti\. On the 
other hand, pure-singlet contributions dominate for x smaller than 0.2. In the limit x — > the 
ratio of twist-4 and twist-2 contributions tends to be a constant. We note that the contribution 
from the quadratic pole turns out to be larger than the contribution from the single one. If 
the pure-singlet quark contribution is replaced by a gluonic one, a very similar shape of twist-4 
correction results, in accordance with the discussion in section 3.2. Figure ^ shows the same 
ratio with a logarithmic scale of the x-axis. 

Fig. 1^ shows hL{x), the renormalon correction to Fl plotted against the coefficient of the 
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Figure 5: Renormalon model prediction for the ratio h2{x) / F2{x, Q'^) for deuteron and proton targets. 
The full line shows the non-singlet and pure-singlet quark contributions. The dashed line corresponds 
to the non-singlet and gluonic contributions. The filled and empty circles denote data points for proton 
and deuteron, respectively [^. The magnitude of the non-singlet contribution has been normalized 
to the large-x data. = 200 MeV, = 4 GeV^ and rif = 4. 



1/Q^ dependent term in the data fit of for a deuteron target. The target mass correction 
taken from [Q, makes up to 15%-20% of hi{x) as compared to the phenomenological fit. The 
pure-singlet part becomes considerable only in the small x region, where it produces a steep 
raise, clearly at variance with the parameterization of Ref. P9 |. 

As expected, in the unpolarized case the pure-singlet contribution is negligible at large x, but 
it dominates when x becomes small. Although, strictly speaking, the absolute normalization is 
unknown, the renormalon model clearly suggests that while the non-singlet twist-4 contribution 
becomes smaller and smaller with x decreasing, the magnitude of the pure-singlet contribution 
raises for x below 0.2, and therefore twist-4 corrections should not be neglected in this region. 

Recall that at large x the UV dominance hypothesis predicts, in accordance with arguments 
based on an analysis of multiparton contribution of hadronic wave functions p3, that power 



corrections to F2{x, Q^) are effectively suppressed by 1/{Q'^{1 — x)) ~ l/W"^, and therefore they 
can be eliminated by a suitable W"^ cut-off. In the small-x domain the situation is different. 
Our calculation predicts that the ratio of twist-4 to twist-2 terms in this region depends rather 
weakly on x, such that the former cannot be minimized by a kinematical cut. 

Here a natural question arises to which extent the renormalon model predictions can be 
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Figure 6: The same as in Fig. but with a logarithmic scale on the x-axis 

trusted in the small x domain. Although it is formulated in the field-theoretical framework, the 
model almost certainly misses some important QCD physics. We note e.g. that ladder corrections 
which are known to infiuence strongly the small x behavior of twist-2 structure functions are 
absent here, as they are formally suppressed by l/n/. Hence, while we expect that general trends 
followed by twist-4 contributions are correctly reproduced, it is possible that details of small x 
behavior of twist-4 corrections are not properly described in this approach. 

Fig. H shows the renormalon model prediction for the twist-4 part of the polarized structure 
function xg^~'^{x,Q'^). The resulting NS correction is negative at small x, changes sign around 
X = 0.02, and reaches a maximum around x ~ 0.4. The pure-singlet contribution is significant 
only in the small x region but, contrary to the unpolarized case, even there it does not dominate 
over the NS part. As in the F2 case the ratio of twist-4 and twist-2 contributions to gi{x,Q'^) 
tends to be a constant for a; — * 0. In order to get a maximal magnitude of the pure-singlet 
correction, in Fig. R| we have counted contributions from single and quadratic poles coherently. 



5 Summary 



The calculation presented in this paper extends previous renormalon model analyses of twist-4 
corrections to DIS [|12], |I3|, |T^, |T^ to the fiavour singlet sector. We have presented a sample of 
predictions for the x-dependence of twist-4 corrections to nucleon structure functions -F2(x, Q^), 
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Figure 7: The shape of the power correction hL{x) to the longitudinal structure function 

Fl{x,Q^) for a deuteron target. The long-dashed line depicts the pure-singlet quark correction, 
with contributions from single and double poles added coherently. The short-dashed line shows the 
contribution of the non-singlet part. The dotted line is the target mass correction, taken from ||50(| . 
The dot-dashed line is the sum of all contributions. The solid line represents hL{x) as obtained 
from the Whitlow-fit to R{x) |^ and the parameterization of from the NMC-analysis p4|| . 
For consistency, the renormalon model prediction shown here has been obtained using the same 
parameterization of twist-2 parton densities. A^jjg = 250 MeV, = 4 GeV^ and Uf = 4. 



Q^), and Fl{x,Q^) in the whole x domain. Our results suggest that in the small-x region 
the twist expansion parameter is not x ■ A^/Q^, as it would follow from a l/VT^-dependence of 
higher twist corrections, but is approximately equal to A^/Q^- In this case higher twists cannot 
be eliminated by a W"^ cut, and should be taken into account in the interpretation of the current 
data. 

One possibility is to fit the data taking into account, in addition to twist-2 contribution, twist- 
4 corrections either in the form derived in the present paper with the non-singlet contribution 
taken from |15|, or as given in Ref. |^ . 

To obtain the renormalon model predictions for twist-4 corrections to nucleon structure 
functions in the flavour- singlet sector one should insert the pure-singlet coefficient functions, 
listed in the appendix below, into the right-hand side of Eq. (^, and add the non-singlet type 
contribution. For each contribution, the overall factor and its mass scale A^ should be fitted 
independently. In particular, mass scales can be different for different channels. In principle one 
should also try a fit in which overall signs differ between contributions arising from single and 
quadratic poles. Because, as discussed in section 3.2, shapes of the resulting quark and gluon 
contributions are rather similar, such a fit could involve, besides standard parameterizations of 
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Figure 8: The shape of the l/Q^ power correction h/^{x) to the polarized structure function 
xgi{x,Q'^) for a deuteron target. The dashed line corresponds to the pure-singlet quark contri- 
bution, the dotted line to the contribution of the non-singlet part. The full line is the sum of both. 

= 231 MeV, = A GeV^ and rif = 4. 

the twist-2 structure functions at some low virtuality, only two additional elements: the quark 
non-singlet term, which can be taken from |T^, and the gluonic term, as given below. Description 
of twist-4 corrections requires then the introduction of 3 new parameters: the normalization of 
the non-singlet quark contribution and the normalization and the mass scale which enters the 
gluonic contribution. In the currently most interesting case of F2{x,Q'^) a further simplification 
is possible: as the flavour non-singlet part is important only in the large-x domain, it can be 
neglected in a fit to the small-x data, which would reduce the number of additional fit parameters 
to 2. 
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6 Appendix 

In this section we quote explicitly our renormalon model results for the pure-singlet coefficient 
functions A^^^ ^'^\z,\og{A'^ /Q'^)) and Af'Q{z,\og{A'^/Q'^)), which control quark and gluon con- 
tributions to (9(1/(5^) corrections to nucleon structure functions in the flavour singlet sector. 
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see Eq. (H). For completeness, we quote also the coefficient functions A^^^ '{z,\og{A'^/Q'^)) and 

A\'^q{z, log(A^/Q^)), arising from the second IR renormalon at s = 2, which can be used to model 
0{1/Q^) corrections. 

To facilitate comparison with the data which treats contributions from single and quadratic 
poles separately, we have splitted each coefficient function into contributions arising from the 
single and quadratic poles, respectively. 



6.1 Quark pure-singlet coefficient functions 



Coefficient function for the quark pure-singlet renormalon contribution to F2: The contribution 
of the quadratic pole is proportional to logi^A"^ z / Q'^) . 



^F—^ 



32, . . 592. 5 ,200 40, , 3 



- + 241og(.)).^ + (— --log(.))) 

32 



\og{A'zlQ' 



z" + (48 \og{z) 



184. 



64\ 

+ (72 + 161og(z))z2__j 



(36) 



<r'^(^,log(AVQ^)) 



Tpn/ 1 
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+ +(— log(^. 75 



, /A 128, 4 ,64 , , , 256, 3 
+ (64 log(2;) + — ) 2;^ + (— log(2;) + — ) z"- 

, 64 , , , 416, 2 512 ,, , 4576 \ 

+ ( log(z) H ) ^ log(;z) 

^5 ' 75' 105 ' 11025/ 



641og(AW)(.^4.^4.^ + 4-^^ + |/ 



(37) 



Coefficient function for the quark pure-singlet renormalon contribution to Fl. The ffist two 
lines contain the contribution from the single pole, the last line contains the contribution of the 
quadratic pole: 



<'^'^(^,log(AVQ^)) 
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(39) 



Coefficient function for the quark pure-singlet renormalon contribution to gi. The contribution 
of the quadratic pole is proportional to log{A'^ z / Q^) . 
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AT/\z,\og{K'/Q')) = 



+(-4-81og(z)>-y 
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6.2 Gluon coefficients functions 



In order to distinguish between contributions arising from single and quadratic poles we have 
split the coefficient functions explicitly into two parts, both proportional to -^wr-- In addition, 
to keep the resulting expressions as compact as possible, we have introduced the notation 



c{z) = \/5cos Q \/71og(2;) 
s(^) = y|sinQ V71og(z) 



(42) 



which appear often in the process of deconvoluting the unpolarized quark-gluon splitting function 
from the pure-singlet coefficient function. 
Gluon coefficient functions for F2: 
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Gluon coefficient functions for Fl: 
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Gluon coefficient functions for gi. 
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